Phase Diagram of the Bose-Hubbard Model with Ta symmetry 
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In this paper we study the quantum phase transition between the Insulating and the globally 
coherent superfluid phases in the Bose-Hubbard model with Ts structure, the "dice lattice". Even 
in the absence of any frustration the superfluid phase is characterized by modulation of the order 
parameter on the different sublattices of the Tz structure. The zero-temperature critical point as a 
function of a magnetic field shows the characteristic "butterfly" form. At fully frustration the super- 
fluid region is strongly suppressed. In addition, due to the existence of the Aharonov-Bohm cages 
at / = 1/2, we find evidence for the existence of an intermediate insulating phase characterized by 
a zero superfluid stiffness but finite compressibility. In this intermediate phase bosons are localized 
due to the external frustration and the topology of the lattice. We name this new phase the 
Aharonov-Bohm (AB) insulator. In the presence of charge frustration the phase diagram acquires 
the typical lobe-structure. The form and hierarchy of the Mott insulating states with fractional 
fillings, is dictated by the particular topology of the Tz lattice. The results presented in this paper 
were obtained by a variety of analytical methods: mean-field and variational techniques to approach 
the phase boundary from the superconducting side, and a strongly coupled expansion appropriate 
for the Mott insulating region. In addition we performed Quantum Monte Carlo simulations of the 
corresponding (2+l)D XY model to corroborate the analytical calculations with a more accurate 
quantitative analysis. We finally discuss experimental realization of the Ts lattice both with optical 
lattices and with Josephson junction arrays. 

PACS numbers: 



I. INTRODUCTION 

The Bose-Hubbard (BH) modeli is a paradigm model 
to study a variety of strongly correlated systems as su- 
perconducting filma^, Josephson junction arrays^ and op- 
tical lattices^. This model predicts the existence of a 
zcro-tcmperature phase transition from an insulating to 
a superfluid state which, by now, has received ample ex- 
perimental confirmation. The BH model is character- 
ized by two energy scales, an on-site repulsion energy 
between the bosons U and an hopping energy t which al- 
lows bosons to delocalize. At zero temperature and in the 
limit U ^ t bosons are localized because of the strong 
local interactions. There is a gap in the spectrum for 
adding (subtracting) a particle, hence the compressibility 
vanishes. This phase is named the Mott insulator. In the 
opposite limit, U <^ t, bosons are delocalized and hence 
are in a superfluid phase. There is a direct transition 
between the Mott insulator and the superfluid state at a 
critical value of the ratio t/U. This Superfluid-Insulator 
(SI) transition has been extensively studied both theoret- 
ically and experimentally and we refer to Refs. 2 3^0 
(and references therein) for an overview of its properties. 

Magnetic frustration can be introduced in the BH- 
model by appropriately changing the phase factors asso- 
ciated to the hopping amplitudes. The presence of frus- 
tration leads to a number of interesting physical effects 
which has been explored both experimentally and theo- 
retically. In Josephson arrays, where this is realized by 
applying an external magnetic field, frustration effects 
have been studied extensively in the past for both clas- 
sical^ and quantum systems^. Very recently a great in- 



terest in studying frustrated optical lattices has emerged 
as well2*2i2iiS. There are already theoretical proposals to 
generate the required phases factors by means of atoms 
with different internal states^ or by applying quadrupolar 
fields^. 

The interest in the properties of dice latticesii has been 
stimulated by the work by Vidal et ali^ on the existence 
of localization, the so called Aharonov-Bohm (AB) cages, 
in fully frustrated dice lattices without any kind of dis- 
order. The existence of these cages is due to destruc- 
tive interference along all paths that particles could walk 
on, when the phase shift around a rhombic plaquette is 
TT. Following the original paper by Vidal et at, several 
experimentaliSiMii^ and theoretical worka^^d''i-^^d^i^°i^-^ 
analyzed the properties of the AB cages. In the case 
of superconducting networks most of the attention has 
been devoted to classical arrays with the exception of 
Refs . IT9II20I where a frustrated quantum quasi-one dimen- 
sional array were studied. 

In quantum arrays (charge) frustration can also be in- 
duced by changing either the chemical potential (in op- 
tical lattices) or by means of a gate voltage (in Joseph- 
son junction arrays). This has the effect of changing the 
electrostatic energy needed to add/remove a boson on a 
given island. The phase diagram present a typical lobe- 
like structursi. Moreover, depending on the range of 
the interaction, it may also induce Wigner-like lattices 
of Cooper pairs commensurate with the underlying lat- 
tice2a. 

The aim of this work is to study the phase diagram of 
a Bose-Hubbard model on a lattice (shown in Fig^). 
We will consider both the cases of electric and magnetic 
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FIG. 1: The Ta lattice: It consist of hubs (with six near- 
est neighbours) connected to rims (three nearest neighbours) . 
The Ta structure is a Bravais lattice with a base inside the 
conventional unitary cell. The lattice vectors are ti and t2. 
The basis is given by the sites A,B,C. Due to the fact that 
these sublattices are not self-connected and have different co- 
ordination numbers, we refer to this structure as tripartite. 
All rhombic plaquettes are identical, although differently ori- 
ented. 



hopping and the Coulomb energy the system undergoes 
two consecutive quantum phase transitions. At the first 
critical point there is a transition from a Mott insulator 
to a Aharonov-Bohm insulator. The stiffness vanishes in 
both phases but the compressibihty is finite only in the 
Aharonov-Bohm insulator. At a second critical point the 
system goes into a superfluid phase. Most of the analysis 
is presented by using approximated analytical methods. 
These results will be checked against Monte Carlo sim- 
ulations that we present in Sec. II V Dl Few details of 
the mapping of the model used in the simulation are re- 
viewed in the Appendix^ The concluding remarks are 
summarized in Sec. 



II. QUANTUM PHASE MODEL ON A % 
ARRAY 



frustration. The location and the properties of the phase 
diagram will be analyzed by a variety of approximate an- 
alytical methods (mean-field, variational Gutzwiller ap- 
proach, strong coupling expansion) and by Monte Carlo 
simulations. The lattice has been experimentally real- 
ized in Josephson arraysMi. In addition we show that it is 
possible to realize it experimentally also with optical lat- 
tices. Although the main properties of the phase diagram 
are common to both experimental realizations, there are 
some differences which are worth to be highlighted. 

The plan of the paper is the following. In the next 
Section we will discuss the appropriate model for both 
the case of a Josephson junction array and optical lat- 
tices (Sec. mi and discuss in some detail how the struc- 
ture can be realized in an optical lattice (Sec. Ill Bp . In 
the same Section we introduce the relevant notation to be 
used in the rest of the paper. A description of the various 
analytical approaches used to obtain the phase diagram 
will be given in Sec. IIIII The zero-temperature phase di- 
agram, in the presence of magnetic and electric frustra- 
tion, will then be described in Sec. IIVI We first discuss 
the unfrustrated case and afterwards we consider the role 
of electric and magnetic frustration respectively. Due to 
the particular topology of the T3 lattice the supercon- 
ducting phase is characterized, even at zero frustration, 
by a modulation of the order parameter on the different 
sublattices (i.e. hubs and rims), which indicates a dif- 
ferent phase localization on islands depending on their 
coordination number. A uniform electrostatic field gives 
rise to a lobe structure in the phase diagram which is 
discussed for the 73 array in Sec. \IVBl The effect of a 
uniform external magnetic field, discussed in Sec. MV C\ 
may induce important qualitative changes in the phase 
diagram in the case of fully frustration. In particular we 
will provide evidences that there is an important signa- 
ture of the Aharonov-Bohm cages in the quantum phase 
diagram. It seems that due to the AB cages a new phase 
intermediate between the Mott insulating and superfluid 
phases should appear. On varying the ratio between the 



Both Josephson arrays and optical lattices are experi- 
mental realizations of the Bose-Hubbard model 

= H-u + Tit 

= ^ X!'^"* ~ no)Uij{ni - no) 

<ij> 

When the mean occupation h on each lattice site is 
large, one is allowed to introduce the phase operator tpi 
by approximating the boson annihilation operator on site 
I by 6i ~ exp[i(^i]. The density rii and phase ipi 
operators are canonically conjugate on each site 

[n„e±*'^^] =±,5,,,e±*'^'. (2) 

In the present work we will focus our attention on the 
quantum rotor version of the model in Eq.JQ) that reads: 

H = ^ ^(rii - no)Z^ij (rij - no) 

- t ^ cos{ipi - ipj - A,^j) . (3) 

The first term on the r.h.s of Eq.Q represents the re- 
pulsion between bosons {lAij depends on the range of 
the interaction and on its detailed form). The second 
term is due to the boson hopping [t = fit is the coupling 
strength) between neighboring sites (indicated with (.) 
in the summation) . The gauge-invariant definition of the 
phase in presence of an external vector potential A and 
fiux-per-plaquette <& ($0 = he/ 2 e is the flux quantum) 
contains the term j = ^ A-dl . All the obscrvablcs 
are function of the frustration parameter dcflned as 
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where the hne integral is performed over the elementary 
plaquette. Due to periodicity of the model it is sufficient 
to consider values of the frustration < / < 1/2. Charge 
frustration is due to a non-integer value hq. As for the 
magnetic frustration also in this case the properties will 
be periodic under the transformation no no+1. Due to 
the additional symmetry uq —hq it is sufficient to con- 
sider value of the charge frustration no in [0, 1/2]. Differ- 
ently from the magnetic frustration the value of n = 1/2 
does not necessarily correspond to fully (charge) frustra- 
tion as this depends on the range of the interaction Ui j. 

The Ta latticsii is represented in Fig^ the lines be- 
tween the sites corresponds to those links where boson 
hopping is allowed. The structure is not itself a Bra- 
vais lattice, but could be considered as a lattice with a 
base inside the conventional unitary cell (see Fig. ^ de- 
fined by the vectors 

ti = (3/2;-V3/2) a 
ta = (3/2; ^-^/3/2) a . 

where a is the lattice constant. The lattice sites of the 
base are at positions 

dA = (0;0)a dB = (0;l)a dc = (0;2)a. 

The reciprocal lattice (ga • tt, = 2TT5a^b) vectors are de- 
fined as 

gi = ^(l/3;-V3/3) 

g2 = ^(l/3;+V3/3) . 

In several situations it turns out to be more convenient 
to label the generic site i by using the position of the cell 
t = niti-fn2t2 {~Ni < ni < Ni) and the position within 
the cell a = A,B,C. In the rest of the paper we either 
use the index i or the pair of labels (t, a). A generic ob- 
servable Wi can be written henceforth as ^^(t). By im- 
posing Born- Von Karman periodic boundary conditions 
its Fourier transform is given by 

W^„(K) ^ , ^ YWa^jt) e-*K-* (5) 

with K = fcigi -I- k2g2 ill the first Brillouin zone. 

It is also useful to introduce a connection matrix T 
whose entries are non-zero only for islands connected by 
the hopping. More precisely Ta^-f{t, t') = 1 if site a of cell 
t is connected by a line (see Fig^J to site 7 of cell t' and 
Ta^-yitjt') = otherwise. The local coordination number 
is thus defined as Za = X^t' -y ^0,7(^7 1'). It is z = 6 for 
the hubs (labelled by A) and z = 3 for the rims (labelled 
by B and C). For later convenience we also define the 
matrix V with elements 

P„,^(t,t') =r,,^(t'-t) e'^".T(t,t') (6) 



which includes the link phase factors which appear if the 
system is frustrated. In the whole paper we fix fc^ = h = 
c 1. 

In the next two subsections we give a brief description 
of the origin and characteristics of the coupling terms in 
the Hamiltonian of Eq. Q both for Josephson and optical 
arrays. In addition we show how to realize optical lattices 
with 73 symmetry. 

A. Josephson junction arrays 

Since the first realization of a Josephson Junctions Ar- 
ray (JJA)S, these systems have been intensively stud- 
ied as ideal model systems to explore a wealth of classi- 
cal phenomena— such as phase transitions, frustration 
effects, classical vortex dynamics and chaos. One of 
the most spectacular result was probably the experimen- 
tal observation^^ of the Berezinskii-Kostcrlitz-Thouless 
transition (BKT)S&. Indeed, well below the BCS transi- 
tion temperature, and in the classical limit, JJAs are 
experimental realization of the XY model. For suffi- 
ciently small (submicron) and highly resistive (normal 
state resistance Rn > Rq = h/Ae^) junctions quantum 
effect start to play an important role. In addition to the 
Josephson energy, which controls the Cooper pair tun- 
nelling between neighboring grains, also the charging en- 
ergy e^/2C (C is the geometrical junction capacitance) 
becomes important. 

Experiments on JJAs are performed well below the 
BCS critical temperature and thus each island is in the 
superconducting state. The only important dynamical 
variable is the phase i^i of the superconducting order 
parameter in each island, canonically conjugated to the 
number of extra Cooper pairs rii present on that island. 
In Eq.Q, the coupling constant t equals the Josephson 
coupling. Hence the second term in Eq.© represent the 
Josephson energy. The first term is due to charging en- 
ergy which can be evaluated by assuming that each island 
has a capacitance to the ground Co and each junction a 
geometrical capacitance C . The electrostatic interaction 
between the Cooper pairs is defined as 

U = 2e^C-^ . (7) 

The capacitance matrix is given by 

= (Co + z,C)4j . (8) 

Since both the connection and the capacitance matrices 
depend only on the distance between the cells (and on 
the base index of both sites), their space dependence can 
be simplified to 

Ca,-yit,t') = Ca,-yiO,t' -t) = Ca,-y{t' -t) (9) 

An estimate of the range o f the electrostatic interaction 
is given bji^i A w y/ C/Cq. The charge frustration uq, 
that we assume to be uniform, can be induced by an 
external (uniform) gate voltage Vq = no/Co- 



4 




FIG. 2: Electrostatic energy (in units of Uo = e^/2Co) re- 
quired to put an extra Cooper pair (for zero external charge) 
on an hub (straight line) and on a rim (dashed) as a function 
of the reduced capacitance C/Co 



Due to the particular structure of the T3 lattice, the 
charging energy of a single (extra) Cooper pair placed 
on a given islands depends on that site being a rim or a 
hub as shown in Fig|21 As a consequence quantum fluc- 
tuations of the phase of the superconducting order pa- 
rameter may be different in the two different cases (rims 
or hubs). We will see in Sec lIV Al that this property is 
responsible for an additional modulation of the order pa- 
rameter in the superconducting phase. 



B. Optical lattices 

Following the work of Jaksch et alX, optical lattices 
have been widely studied as concrete realization of the 
Bose-Hubbard model that is, as we saw, directly re- 
lated to the quantum phase model studied in this paper. 
The experimental test of the SI transition^ has finally 
opened the way to study strongly correlated phenom- 
ena in trapped cold atomic gases. Very recently, several 
works addressed the possibility to induce frustration in 
optical Iattice3£i2i2ii£. It is therefore appealing to test the 
properties of the lattice also with optical lattices once 
it is known how to create 7^ lattices by optical means. 

Here we propose an optical realization of a ^3 struc- 
ture by means of three counter-propagating pairs of laser 
beams. These beams divide the plane in six sectors of 
width 60 (see the inset of FiglSJ and are linearly po- 
larized such to have the electrical field in the xy plane. 
They are identical in form, apart from rotations, and 
have wavelength equal to A = 3/2 a (a is the lattice con- 
stant. Given a polarization of a pair of lasers on the 
y-axis El = {Q,Ey) the other two pairs are obtained by 
rotating Ei of 120 and —120 around the z-scxis. The 
square modulus of the total field gives rise to the desired 




FIG. 3: Optical potential with Tg symmetry generated by 
three counter-propagating laser beams. The inset shows the 
bidimensional contour plot while in the figure the details of 
the profile along a line connecting three sites (placed at posi- 
tions X = 1, X = Q and x = —1) is shown. The sites x = —1, 1 
are rims while the site at a; = is a hub. Also here, as in the 
case of JJAs, the different form of the potential implies that 
the on-site energy i7o is different for hubs and rims. 

optical potential as it is shown in Fig|31 

The form of the potential landscape also in this case 
imposes that the on-site repulsion may be different for 
hubs and rims. It is however diagonal 

U ^Urlr + Uh lh ■ (10) 

The subscript /i, r denotes the respectively the hub and 
rim sites and Th^r are the projectors on the corresponding 
sublattices. In Eq.©, now, the coupling t describes the 
hopping amplitudes for bosons, no is proportional to the 
chemical potential, Aij is the effective "magnetic frustra- 
tion" that in this case may have several different origin 
depending on the scheme used. For simplicity we will 
always refer to A as to the vector potential and we will 
use the magnetic picture also for optical lattices. 



III. ANALYTIC APPROACHES 

The SI transition has been studied by a variety of 
methods; here we apply several of them to understand 
the peculiarities that emerge in the phase diagram due 
to the 79 lattice structure. The results that derive from 
these approaches will be presented in the next section. 

The location of the critical point depends on the exact 
form and the range of hiij. This issue is particulary in- 
teresting when discussing the role of electric frustration. 
In the paper we address the dependence of the phase 
boundary on the range of the interaction in the mean- 
field approximation. The variational Gutzwiller ansatz 
and the strong coupling expansion will be analyzed only 
for the on-site case of Ea. (|10|l . In the case of magnetic 
frustration the form of lAi^j leads only to quantitative 
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changes so, also in this case, we discuss only the on-site 
case. 



A. Mean field approach 

The simplest possible approach to study the SI phase 
boundary consists in the evaluation of the superconduct- 
ing order parameter, defined as 

by means of a mean-field approximation. By neglect- 
ing terms quadratic in the fiuctuations around the mean 
field value, the hopping part of the Hamiltonian can be 
approximated as 

The order parameter is then determined via the self- 
consistency condition 

Tr |e*^>(-') e-/3«^ T, exp (^j^ 7^^'"^'(t)) } 

"^'^^ ^ " Tr{e-^'«-T.exp(/>(™^)(r))} 

(11) 

In the previous equation, Tr is the time-ordering in imag- 
inary time r and (3 — 1/T . The r dependence of 
the operators is given in the interaction representation 
W{t) = e"^^^ W e~'^ . For simplicity we already as- 
sumed the order parameter independent on the imagi- 
nary time. One can indeed verify that this is the case in 
the mean-field approximation. Close to the phase bound- 
ary the r.h.s. of Eq.^]can be expanded in powers of the 
order parameter and the phase boundary is readily de- 
termined. 

A central quantity in the determining the transition is 
the phase-phase correlator 

G^jM = (Tre*"^'(^)e-"^^("V (12) 

where the average is performed with the charging part of 
the Hamiltonian only. Charge conservation imposes that 
the indexes i,j are equal. The Matsubara transform at 
T = of the correlator reads 

/OO 1 
G,,(r)e'- = 5] — (13) 

where /S.Ea^± are the excitation energies (for zero Joseph- 
son tunnelling) to create a particle (-I-) or a hole (-) on a 
site of the sublattice a where i lies. 

In the case of the Ta lattice considered here even at 
zero magnetic field the order parameter is not uniform. 
The tripartite nature of the lattice results in a vectorial 
mean field if) with one component for each sublattice. In 
the general case the linearized form of Ea. ljlll) can be 
rewritten as 

V'o(t) = ^^5]G„,„(0)P„,^(t,t')V7(t') (14) 

7 t' 



that, due to the topology of the lattice is equivalent to 

7 t' 

The phase transition is identified with a non-trivial 
solution to this secular problem, i.e. one should deter- 
mine TTmax, the largest eigenvalue of P. This requirement 
translates in the following equation for the critical point 



ter = 2 = (15) 

In deriving the previous equation we used the fact that 
the sites B and C in the elementary cell (see Fignj have 
the same coordination number and therefore the phase- 
phase correlator is the same. In addition to the evalua- 
tion of the Matsubara transform at zero frequency of the 
phase correlator, one has to determine the eigenvalues of 
the gauge-link matrix V. With a proper gauge choice it 
is possible to reduce this matrix to a block diagonal form. 
For rational values of the frustration, / = p/q, by choos- 
ing A = {x — \/3 y) f y, the magnetic phase 
factors Aij{t,t') (shown in Fig. 0)) have a periodicity 
of r X 1 elementary cells with r = LCM{q,3)/3. This 
implies that in the Fourier space (see Eq.©) the com- 
ponent /c2 is conserved and that ki is coupled only with 
the wavevectors k["^^ = ki + 2'Km/r (m = 0, ...r — 1). 
The determination of i^max is therefore reduced to the 




FIG. 4: Magnetic phase pattern with the gauge choice 
A = - V3y)^/y 



diagonalization of a 3r x 3r matrix (Pq, .^(fci) is r x r) 

/ PAAkl) PAc(fcl) \ 

P(fci,fe) =4.,o PaA^i) 

\p\Aki) ; 

(16) 

with (fci,fc2) belonging to the reduced Brillouin zone 
B,. = {0 < fc, < 27r/r}. 

The matrix V has r zero eigenvalues, and r pairs of 
eigenvalues equal in absolute value given by the reduced 
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secular equation 

Pl,2ikl)Pl2{kl) + Fl,3(fci)P|,3(fci) 



Vl = TT Vl 



This simplification allows us to deal with r x r matrices 
instead of q x q. 

The inclusion of a finite range interaction, important 
only for Josephson arrays, leads to a richer lobe structure 
in presence of electrostatic frustration. The calculation 
of the lobes will be done within the mean field theory 
only. 



B. Gutzwiller variational approach 

A different approach, still mean-field in spirit, that al- 
lows to study the properties of the superconducting phase 
is the Gutzwiller variational ansatz adapted to the Bose- 
Hubbard model by Rokhsar and KotliaB^S The idea is 
to construct a variational wave-function for the ground 
state starting from the knowledge of the wave-function in 
the absence of the interaction term Tiu in the Hamilto- 
nian. In this case, and in absence of magnetic frustration, 
the ground state has all the phases aligned along a fixed 
direction 9. In the boson number representation it reads 



\GS) 



\{n^}) 



(17) 



A finite charging energy, tends to suppress the compo- 
nents of the state with large charge states, a variational 
state can then be constructed through the ansatz 



1^5) = c„i,...,„„|{n,;}) 



where 



(18) 



(19) 



In Ea. (|19|l Nqs is a normalization factor and Ki and 
Hi are variational parameter to be determined by mini- 
mizing the ground state energy. The Mott insulator is 
characterized hy K = oo, i.e. by perfect localization of 
the charges, K = is the limit of zero charging, a fi- 
nite value of K describes a superfluid phase where the 
phase coherence has been established albeit suppressed 
by quantum fluctuations. 



C. Strong coupling perturbation theory 



Both methods illustrated in Sections llII Al and llll Bl are 

based on the analysis of the superconducting phase and 
on the determination of the phase boundary as the loca- 
tion of points where the superfluid order parameter van- 
ishes. A complementary approach, which analyzes the 



phase boundary from the insulating side, was developed 
by Freericks and Monien^S. The method was applied to 
the case of square and triangular lattices in Ref. [s^ for 
the Bose- Hubbard model and in Ref. |^ for the quantum 
rotor model. In this section we describe how to adapt 
the method to the Ta lattice. We will present the re- 
sults of this analysis, particularly important for the fully 
frustrated case, in Sec. IIV CI 

In the insulating phase the first excited state is sepa- 
rated by the ground state by a (Mott) gap. In the limit of 
vanishing hopping the gap is determined by the charging 
energy needed to place/remove an extra boson at a given 
lattice site. The presence of a finite hopping renormalizes 
the Mott gap which, at a given critical value, vanishes. 
The system becomes compressible, and the bosons, since 
are delocalized, will condense onto a superfluid phase. It 
is worth to emphasize that the identification of the SI 
boundary with the point at which the gap vanishes is 
possible as the bosons delocalize once the energy gap is 
zero. As we will see, in the case of lattice the situation 
becomes more complex. In the presence of external mag- 
netic frustration it may happen that though the Mott gap 
is zero, the states are localized and therefore the charges 
cannot Bose condense. In this cases between the Mott 
and superconducting region an additional compressible 
region (with zero superfluid stiffness) may appear. In 
order to keep the expressions as simple as possible we 
consider only the case of on-site interaction, though we 
allow a different U for hubs and rims as in Ea. (|ll)|l . The 
possible existence of such a phase, however, does not de- 
pend on the exact form of Ui_j. The strong coupling 
expansion is particularly useful for 7^ lattice as it may 
help in detecting, if it does exist, the intermediate phase. 

In the strong-coupling approach of Freericks and 
Monien the task is to evaluate, by a perturbation expan- 
sion in t/U, the energy of the ground and the first ex- 
cited state in order to determine the point where the gap 
vanishes. We denote the ground and first excited levels 
by El''j and Eff^ respectively. The choice of the starting 
point for the perturbation expansion is guided by the na- 
ture of the low-lying states of the charging Hamiltonian. 
When no < 1/2 (and in zero-th order in t/U) the ground 
state of the electrostatic Hamiltonian is (n^ = Vi) and 
first excited level is given by a single extra charge local- 
ized on a site. Levels corresponding to charging a hub and 
a rim are nearly degenerate (i.e. [Ur—Uh) / {Ur+Uh) ^ 1, 
with the hub being lower in energy). As the strength of 
the hopping is increased, the insulating gap decreases. 
We would like to stress, and this is an important differ- 
ence emerging from the Ta topology, i.e. the location of 
the extra charge (on a hub or a rim) requires a different 
energy. This in turn has important consequences in the 
structure of the perturbation expansion. 

Up to the second-order in the tunnelling, the ground 
state energy at tiq = is given by 



p9s 



2 • 2N 



iUh + Ur)/2 4 



(20) 
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FIG. 5: Intermediate charge states involved in the definition 
of Eg. 1251 . In the upper panel the contributions to the di- 
agonal part are shown while in the lower panel there are the 
contributions to the off-diagonal part. The processes repre- 
sented here are those contributing to the second order in the 
hopping amplitude. The black/white circles represent one ex- 
tra -I-/- Cooper pair on a given site. 



where N is the number of sites and 2N the number of 
hub-rim Unks in the lattice. Note that the first-order 
correction vanishes because the tunnelling term does not 
conserve local number of particles. 

Due to nearly degeneracy of the excited levels, one is 
not allowed to perturb each of them independently but 
has to diagonalize the zeroth and the first order terms 
simultaneously. One has to diagonalize the following ma- 
trix: 



(21) 



This task can be reduced to the diagonalization of a 
3r(/) X 3r(/) matrix with a proper choice of the gauge 
(see Section UlI A|) . 

For example, the (degenerate) lowest eigenvalue at / = 
1/2 is 



Q 



(1) 



Uh + Ur 



/=l/2 



(22) 



which reduces to U/2 — t\/6/2 in the case of perfectly 
degenerate charging energy. It must be stressed that all 
the energy bands are flat, independently of the values of 
the charging energies (it depends only on the peculiar V 
structure). 



The second order perturbation term should be calcu- 
lated on the lowest energy manifold: moreover only ma- 
trix elements between states of the same manifold are 
allowed. Nonetheless, it is simpler to write the differ- 
ent contributions in the usual basis of hub and rims (see 
Fig-©)- The first excited state, to second order in tun- 
nelling is given by 



E 



(1) _ 



l(2) 



(23) 



where Q'^' is the second order matrix and can be split 
into separate sub-matrices on different sub-lattices, i.e. 



Q 



(2) 



(24) 



Such a decomposition is possible because after two tun- 
nelling events the boson come back to the initial sublat- 
tice. 



2(2);. 



{Uh -Ur)/2 

+ {2-2N~2zh) 
, - Zhlh 



Zh 



{Uh-{^Uh + Ur))/2 



{Uh - {2Uh + Ur))/2 

{Uh - Ur)/2 ^ {Uh - i2Uh + Ur)) 12 



(25) 



(Q(2)r is defined in a similar way) where X/j are the 
projectors on the hub and rim sublattices. After some 
algebra and by changing basis to the one composed by 
the eigenvectors of Eq. , one gets the first excited en- 
ergy level. The task is now to determine the location 
of points at which the gap, given by the difference of 
Eg. lP^ and Ea. (pn|) . vanishes. It is worth to stress that 
the thermodynamically divergent contributions wash out 
exactly their analogous in the ground state expression of 
EqEOI 

We discuss the results deriving from this approach in 
the next Section where we analyze the phase diagram. 



IV. PHASE DIAGRAM 

In order to keep the presentation as clear as possible 
we first discuss the main features of the phase diagram 
by means of the analytical approaches introduced before. 
We will then corroborate these results in a separate sec- 
tion by means of the Monte Carlo simulations. 

The value of the critical Josephson coupling as a func- 
tion of the range of the electrostatic interaction, in the 
absence of both electric and magnetic frustration is dis- 
cussed first. The effect of frustration, either electric or 
magnetic will then be discussed in two separate sections. 
In the case of electrical frustration the topology of a Ta 
lattice gives rise to a rather rich lobe structure, the over- 
all picture is nevertheless very similar to the one encoun- 
tered in the square lattice. Much more interesting, as 
one would suspect, is the behaviour of the system as a 
function of the magnetic frustration. The location of the 
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phase boundary shows the characteristic butterfly shape 
with an upturn at fully frustration typical of the 7^. In 
addition, at / = 1/2, a very interesting point which 
emerges from our analysis is the possibility of an interme- 
diate phase, the Aharonov-Bohm insulating phase, sepa- 
rating the Mott insulator from the superfluid. 

A. Zero magnetic & electric frustration 

A first estimate for the location of the phase bound- 
ary can be obtained by means of the mean-field ap- 
proach described in Sec. IIII Al The results coincide 
with the first-order perturbative calculation introduced 
in Sec. IIII Cl and with the Gutzwillcr variational approach 
of Sec. IIII Bl In absence of frustration the K = mode 
corresponds to the maximum eigenvalue of the matrix V 
i'n'max = \/l8) and the transition point is given by 

t,r = ^ ^JUAA0)UB.B{0) (26) 

In the limit of on-site uniform {Ur — Uh ~ SUq) the SI 
transition occurs at the value tcr/Uo = 2-\/2/3 ~ 0.943 
very close to the mean field value for a square lattice 
tcr/Uo = 1 (in both lattices the average value of nearest 
neighbours is 4). In the case of a Josephson array the 
transition point depends on the range of the interaction, 
see Eq.®. In the (more realistic) case of a finite junc- 
tion capacitance an analytic form is not available and the 
numerical phase boundary is shown in Fig. as a func- 
tion of the ratio C/Cq. In the case of optical lattices, 




FIG. 6: Josephson arrays; dependence of the critical point 
at / = on the range of the Coulomb repulsion determined 
by the ratio C/Cq. 

see Ea. (|l()|) . the repulsion is on-site. There is still a weak 
dependence of the transition on the difference Ur — Uh- 
As it is shown in Fig. [7| this dependence is not partic- 
ularly interesting and in the Monte Carlo simulation we 
will ignore it. 

As already mentioned, a characteristic feature that 
emerges in T3 lattices, even in the absence of magnetic 




FIG. 7: Optical lattices: Dependence of the transition point 
on the difference repulsion in the hubs and the rims 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 8: Modulation of the order parameter for zero frus- 
tration, V'hub (straight) is always higher than ^prim (dashed). 
The curves are obtained by means of the Gutzwiller varia- 
tional approach. 

frustration, is that the superfluid order parameter is not 
homogeneous. This can be already seen from the eigen- 
vector corresponding to the solution of Ea. (|26|l . Near the 
transition point the ratio between the order parameter 
value on hubs and rims is constant and is related to the 
ratio of the on-site repulsions \^ph/'|pr\ — \/ziJJrJzrUh- 
Phase localization is more robust on hubs (zh = 6) than 
on rims (z^ = 3) because of the larger number of nearest 
neighbours. In order to better understand the modula- 
tion of the order parameter we analyzed the properties 
of the superconducting phase using the variational ap- 
proach exposed in Sec. IIII Bl (which allows us to study 
the behaviour of tJj also far from the transition). As it can 
be clearly seen from FiglHl quantum fluctuations have a 
stronger effect on the rims than hubs due to the different 
coordination number of the two sublattices. Note that 
this is a pure quantum mechanical effect, in the classical 
regime all phases are well defined and ^phub = "iprim — 1- 
The transition point (as it was implicit in the previous 
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3x3 super-cell turns out to be sufficient. Given a certain 
rational filling p/ q, the corresponding charging energy is 
given by 



FIG. 9: Ground state configurations of the charges (i.e. at 
t = 0) for filling 1/9,2/9,1/3,4/9. The different ground 
states occurs on increasing the value of the external charge 
no The black circles denote those sites that are occupied by 
one Cooper pair. The ground state configurations are respon- 
sible for the behaviour of the phase correlator and hence of 
the lobe-like structure, Fig llOl of the phase diagram. 



discussion) is the same for both sublattices: there is no 
possibility to establish phase coherence between rims if 
the hub- network was already disordered (and viceversa). 



B. Electric frustration 

When an external uniform charge frustration is 
present, the array cannot minimize the energy on each 
site separately, hence frustration arises. The behaviour 
of the transition point as a function of the offset charge 
shows a typical lobe-structurei*S^. At the mean-field 
level all the information to obtain the dependence of 
the phase boundary on the chemical potential (gate po- 
tential for Josephson arrays) is contained in the zero- 
frequency transform of the Green functions G in Ea. (|15|l . 
The calculation of the phase-phase correlators, defined by 
Ea. (|12|l . is determined, at T = 0, once the ground and 
the first excited states of Tiu is known. As all the ob- 
servables are periodic of period one in the offset charge 
no and are symmetric around no = 0, the analysis can 
be restricted to the interval [0, 1/2]. Ground state charge 
configuration in the case of some values of the electric 
frustration are shown in Fig|3 

The phase diagram in the presence of charge frustra- 
tion has a lobe structurei in which, progressively on in- 
creasing the external charge, the filling factor increases 
as well. In the case of finite range charging interaction 
also Mott lobes with fractional fillings appear—. An ana- 
lytical determination of the ground state of the charging 
Hamiltonian for generic values of the external charge is 
not available. We considered rational fillings of the whole 
lattice as made up of periodic repetitions of a partially 
filled super-cell of size comparable with the range of the 
interaction Uij and then constructed a Wigner crystal for 
the Cooper pairs with this periodicity. For C /Cq < 1 a 



= 3iV 



Co 



Uq - 2-no 

q 



Co 
N 



E 



where A'^ is the number of cells in the system and {n^} is 
the particular realization of the filling. 

This defines a set of parabolas which allow to deter- 
mine the sequence of ground states. The variation of the 
ground state configurations as a function of gate charge 
gives to the phase boundary a characteristic structure 
made of lobes, as shown in Fig. ^] The longer is this 
range of the electrostatic interaction the richer is the lobe 
structure. 

As can be seen in Fig. llOl when the interaction is purely 
on-site there is only one lobe that closes at half filling 
when the degeneracy between the empty ground state 
and the extra-charged one leads to superconductivity for 
arbitrarily small t. As soon as the range becomes finite, 
other fillings come into play. An interesting feature typ- 
ical of the Ts lattice is that at no = 1/2 the half filled 
state is not the ground state (see Fig. I10|l . 

Finally, we recall that the presence of the offset breaks 
the particle-hole symmetry and thus the universality 
class of the phase transition changai. This can be seen 
from the expansion at small uj of the correlator (Eq. I13|) 
that enters the quadratic term of the Wilson- Ginzburg- 
Landau functional. With no also terms linear in uj en- 
ter the expansion and the dynamical exponent z changes 
from 1 to 2. 



C. Magnetic frustration and Aharonov-Bohm 
insulating phase 

The outgrowing interest in 7^ lattices is especially due 
to their behaviour in the presence of an externally ap- 
plied magnetic field. The presence of a magnetic field 
defines a new length scale, the magnetic length. The 
competition between this length and the lattice period- 
icity generates interesting phenomena such as the rising 
of a fractal spectrum a la Hofstadter. In lattices per- 
haps the most striking feature is the complete localiza- 
tion in a fully frustrating field (/ = 1/2). This is due 
to destructive interference along all paths that particles 
could walk on, when the phase shift around a rhombic 
plaquette is tt (see Fig. lll|l . Is there any signature of this 
localization (originally predicted for tight-binding mod- 
els) in the quantum phases transition between the Mott 
and the superconducting phases? This is what we want 
to investigate in this section. 

In order to determine the phase boundary at T = we 
can follow either the mean field approach of Sec. lIII Al or 
the perturbative theory presented in Sec. IIII CI Wc re- 
mind that while the first approach signals the disappear- 
ance of the superfluid phase, the perturbation expansion 
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FIG. 10: Lobe structures at different values of the capacity, i.e. electrostatic range (respectively C — 0, 10~^, 10~^, 1). The 
dashed lines point out the discrete filling of the ground state. Pictures on the right are magnifications of the highlighted areas 
in the left ones. 




FIG. 11: Aharonov-Bohm cages. Particles that starts on 
white sites can't go further than black sites, due to destructive 
interference. In fact, f — 1/2 means tt phase shift around a 
plaquette. In square lattices this could not happen because 
of the escape opportunity given by straight lines. 



indicates where the Mott phase ends. The results of both 
approaches are shown in Fig. E| Commensurate effects 
are visible in the phase boundary of Fig. ^] at ratio- 
nal fractions f = p/q oi the frustration. The results 
presented are quite generic. We decided to show, as a 
representative example, the results for a JJ array with 
capacitance ratio C /Cq ~ 1 and an optical lattice with 
Ur — Uh = 0.5Uh- The peak at / = 1/2, characteristic 
of the 79 lattice is due to the presence of the Aharonov- 
Bohm cages. 

Although there is a difference between the mean-field 
and the strong coupling calculation, they both confirm 



the same behaviour. A very interesting point however 
emerges at half-filling. It is worth to stress again that 
while the mean-field shows the disappearance of the su- 
perconducting phase, the strong coupling expansion in- 
dicates where the Mott gap vanishes and hence charges 
can condense. The vanishing of the gap can be associ- 
ated to boson condensation only if bosons are dclocalizcd. 
This is the case for the whole range of frustrations except 
at / = 1/2. In the fully frustrated case the excitation 
gap vanishes but the excited state (the extra boson on 
a hub) still remains localized due to the existence of the 
Aharonov-Bohm cages. This may lead to the conclusion 
that at fully frustration there is an intermediate phase 
where the system is compressible (the Mott gap has been 
reduced to zero) with zero superfluid density (the bosons 
are localized in the Aharonov-Bohm cages). 

At this level of approximation there is no way to ex- 
plore further this scenario. In order to assess the exis- 
tence of the intermediate phase a more accurate location 
of the phase boundaries is necessary. We will discuss 
the possible existence of the Aharonov-Bohm insulator 
by means of Monte Carlo simulations in the next section. 



D. MonteCarlo methods 

The simulations are performed on an effective classical 
model obtained after mapping the model of Eq. (|3J) onto 
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FIG. 12: Phase boundary in presence of a magnetic field in Tij: 
straiglit line is the perturbative result, mean field is dashed. 
Upper: JJAs with C/Co = 1; lower: optical lattices with 
Ur — Uh = 0.5Uh- Note the highly pronounced peak at / = 1/2 
in contrast to the square lattice case. 



a (2 + 1) XY model. Our main interest in performing the 
Monte Carlo simulation is to look for signatures of the 
Aharonov-Bohm insulator. As its existence should not 
depend on the exact form of the repulsion Uij we chose 
the simplest possible case in which the repulsion is on-site 
and Uh ^ Ur- The details of the mapping are described 
in Refs. I32l33l and are briefly reviewed in Sec. ^ The 
effective action S (at zero charge frustration) describing 
the equivalent classical model is 

S = K ^ [1 - cos {ipi^k ~ fj.k - Aij)] 

+ K i^~cos{ip,^k~V^M')) ■ (27) 

i,{k, k') 

where the coupling K is y^t/U. The index k labels the ex- 
tra (imaginary time) direction which takes into account 
the quantum fluctuations. The simulations where per- 
formed on L X L X Lr lattice with periodic boundary 
conditions. The two correlation lengths (along the space 
and time directions) are related by the dynamical expo- 
nent z through the relation ^ C^- For zero magnetic 



frustration, because of the particle-hole symmetry (we 
consider only the case no = 0) holds z = 1. As we will 
see this seems not to be the case at fully frustration be- 
cause of the presence of the Aharonov-Bohm cages. 

The evaluation of the various quantities have been ob- 
tained averaging up to 3 x 10^ Monte Carlo configurations 
for each one of the 10^ initial conditions, by using a stan- 
dard Metropolis algorithm. Typically the first 10^ were 
used for thermalization. The largest lattice studied was 
24 X 16 X 24 at fully frustration and 48 x 48 x 48 at / = 0. 
This difference is due to the much larger statistics which 
is needed to obtain sufficiently reliable data. While in the 
unfrustrated case we took a cube of length L in the fully 
frustrated case it turned out to be more convenient to 
consider (but will discuss other lattice shapes) an aspect 
ratio of 2/3. With this choice the equilibration was sim- 
pler probably due to a different proliferation of domain 
wallsi&ii. 

In order to characterize the phase diagram we stud- 
ied the superfluid stiffness and the compressibility of the 
Bose-Hubbard model on a Ta lattice. The compress- 
ibility, K, is defined by k = d^T/dV^ where T is the 
free energy of the system and Vq the chemical potential 
for the bosons. By employing the Josephson relation in 
imaginary time, see Refi^, the compressibility can be 
expressed as the response of the system to a twist in 
imaginary time, ipi^k ^i,k + Or k, i.e. 



del 



(28) 



The superfluid stiffness is associated to the free energy 
cost to impose a phase twist in a direction e, i.e. (fi 
ifi + ■ Vi , through the array 



del 



(29) 



B = 



1. f^O 

In the case of unfrustrated system we expect that the 
transition belongs to the 3D — XY universality class. 
Close to the quantum critical point k ^ £,~^, the corre- 
sponding finite size scaling expression for the compress- 
ibility reads 

n = L-('-^^n[L'f^^^,^) (30) 



An analogous expression holds for the finite size-scaling 
behaviour of the stiffness 



(31) 



The expected exponent is i/ = 2/3 as it is known from 
the properties of the three-dimensional XY model. 
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FIG. 13: / = 0: a)Scaling (main) and data collapse (inset) of 
the compressibility for the unfrustrated case. b)The same as 
in a) for the superfluid stiffness. All the systems have aspect 
ratio Lt = Ly = with L-r = & (circles), 12 (squares), 18 
(diamonds), 24 (triangles up), 30 (triangles down) 



The results of the simulations for the compressibility 
and for the stiffness are reported in Fig^l Finite size 
scaling shows that the SI transition occurs at 

= 0.435 ± 0.0025 . (32) 

As expected the unfrustrated case follows remarkably 
well the standard picture of the Superfluid-Mott Insu- 
lator quantum phase transition. In the absence of the 
magnetic field the system defined by Ea. l|77|l is isotropic 
in space-time and therefore the stiffness and the com- 
pressibility have the same scaling and critical point. 



a- S = \ 

The situation changes dramatically in the fully frus- 
trated system. In this case an anisotropy in space and 
time directions arises because of the presence of the ap- 
plied magnetic field which frustrates the bonds in the 
space directions (see the r.h.s of Ea. (|27|l 'l. This field 



induced anisotropy is responsible for the different be- 
haviour of the system to a twist in the time (compress- 
ibility) or space (stiffness) components. 

As already observed in the classical casaii, the Monte 
Carlo dynamics of frustrated systems becomes very 
slow. This seems to be associated to the presence of 
zero-euCTgy domain walls first discussed by Korshunov 
in Rcf. [l^ This issue is particulary delicate for the su- 
perfluid stiffness. In this case the longest simulations 
had to be performed. Moreover in order to alleviate this 
problem we always started the run deep in the super- 
fluid state and progressively increased the value of the 
Hubbard repulsion U . Also the choice of the lattice di- 
mensions turned out to be important. We made the sim- 
ulations on 12 X 8 X 12, 18 x 12 x 18, and 24 x 16 x 24 
systems and found out that by choosing this aspect ratio 
along the x and y directions thermalization was consid- 
erably improved. 

The results of the simulations are reported in Figll4l 
for the compressibility and for the stiffness. As it appears 
from the raw data of the figure it seems that the points 
at which the compressibility and the stiffness go to zero 
are different. An appropriate way to extract the critical 
point(s) should be by means of finite size scaling. 

As a first attempt we assumed that the transition is in 
the same universality class as for the unfrustrated case 
and we scaled the data as in Fig^l Although the scal- 
ing hinted at the existence of two different critical points 
for the Mott to Aharonov-Bohm insulator and for the 
Aharonov-Bohm insulator to superfluid transitions re- 
spectively, the quality of the scaling points was poor. In 
our opinion this observation may suggest that the scaling 
exponents for the fully frustrated case are different as the 
one for the direct Mott Insulator to Superfluid phase 
transition at / ^ 1/2. In order to extract more tight 
bounds on the existence of this phase we analyzed the 
size dependence of the observables without any explicit 
hypothesis on the scaling exponent (which we actually 
do not know). The results are presented in Fig llSI The 
data of Fig^lseem to indicate that there is a window 

0.65 < K-'^ < 0.7 

where the system is compressible but not superfluid! This 
new phase, the Aharonov-Bohm insulator, is the result 
of the subtle interplay of the T3 lattice structure and 
the frustration induced by the external magnetic field. 
Our simulations cannot firmly determine the existence of 
two separate critical points since we were not able to im- 
prove their accuracy and study larger lattices. However 
we think that, by combining both the analytical results 
and the Monte Carlo data we have a possible scenario for 
the phase diagram of the frustrated BH model on a Ta 
lattice. 

Further evidence of the existence of the AB cages can 
be obtained by analyzing the anisotropy in space and 
time directions of the phase correlations. For this pur- 
pose we considered the compressibility as a function of 
L and separately. The idea is that because of the 
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FIG. 14: / = l/2:Compressibility (upper panel) and stiffness 
(lower panel) assuming the aspect ratio Lx Lx2L/3. Different 
symbols corresponds to L = 12 (circles), 18 (squares), 24 
(triangles up), and 30 (triangles down). 



AB cages the correlations are short-ranged in the space 
directions (bosons are localized) while there are longer 
ranged correlations in the time direction. Indeed the de- 
pendence of the compressibility on the system dimen- 
sions is strong when one changes Lr while it is rather 
weak when the space dimensions are varied as shown in 
FigEl This hints at the fact that the Aharonov-Bohm 
phase is a phase in which the gap has been suppressed 
(correlation in the time dimension) but where the bosons 
are localized (short-range correlations in space). 

The Monte Carlo simulations just discussed provide ev- 
idence for the existence of a new phase between the Mott 
insulator and superfluid. Due to the finite size of the 
system considered and to the (present) lack of a scaling 
theory of the two transitions, wc cannot rule out other 
possible interpretations of the observed behaviour of the 
Monte Carlo data. A possible scenario which is com- 
patible with the simulations (but not with the result of 
the perturbation expansion^) is that a single thermo- 
dynamic transition is present in the 2-1-1 dimensional 
system but the phase coherence is established in a two 
step process. First the system becomes (quasi) ordered 
along the time direction, then, upon increasing the hop- 




FIG. 15: / = 1/2: Compressibility (upper panel) and the stiff- 
ness (lower panel) as a function of the Lt size of the system for 
different values of K. Data corresponds to 1/K = 0.6(circles), 
0.625(squares), 0.65(diamonds), 0.675(triangles up), 0.7(tr. 
left), 0.725(tr. down), 0.75(tr. right). 
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FIG. 16: Compressibility as a function oi K — \J t/U for 
different values of the system sizes. Different symbols cor- 
responds to different lengths L-r in the time dimension: 6 
(circles), 12 (squares), and 30 (diamonds). Different fillings 
are different spatial sizes Lx x Ly-. 6x4 (black), 12 x 8 (gray), 
and 18 x 12 (white) . The compressibility depends strongly on 
Lr but very weakly on Lx x Ly. 
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ping the residual interaction between these "quasi-one- 
dimensional" coherent tubes go into a three-dimensional 
coherent state driven by the residual coupling between 
the tubes. In more physical terms the "tubes" repre- 
sent the boson localized in the AB cages and the residual 
hopping is responsible for the transition to the superfluid 
state. This means that the intermediate state that we ob- 
serve is due to a one- to three-dimensional crossover that 
takes place at intermediate couplings. 



V. CONCLUSIONS 

In this work we exploited several methods, both an- 
alytic and numerical, in order to determine the phase 
diagram of a Bose-Hubbard model on a Tij lattice. Dif- 
ferently from previous studies on 7^ networks we ana- 
lyzed the situation where the repulsion between bosons 
(or Cooper pairs for Josephson arrays) becomes compara- 
ble with the tunnelling amplitude (Josephson coupling in 
J JAs) leading to a quantum phase transition in the phase 
diagram. Up to now the attention on experimental imple- 
mentations has been confined to Josephson networks. As 
discussed in Sec. Ill Bl the Ta lattice can also be realized in 
optical lattices. The possibility to experimentally study 
frustrated 7^ optical lattices open the very interesting 
possibility to observe subtle interference phenomena as- 
sociated to Aharonov-Bohm cages also with cold atoms. 
Having in mind both the realization in Josephson and 
optical arrays, we studied a variety of different situations 
determined by the range of the boson repulsion including 
both electric and magnetic frustration. Although in the 
whole paper we concentrated on the T = case, in this 
discussion we will also comment on the finite temperature 
phase diagram. 

The peculiarity of the lattice symmetry already 
emerges for the unfrustrated case. The superfluid phase 
is not uniform but it has a modulation related to the pres- 
ence of hubs and rims with different coordination num- 
ber. As a function of the chemical potential (gate charge) 
the transition has a quite rich structure due to the dif- 
ferent boson super-lattices which appear as the ground 
state. 

As a function of the magnetic field the SI transition has 
the characteristic butterfly form. In the fully frustrated 
case, however, the change is radical and we find indica- 
tions that the presence of the Aharonov-Bohm cages can 
lead to the appearance of a new phase, the Aharonov- 
Bohm insulator. This phase should be characterized by 
a finite compressibility and zero superfluid stiffness. A 
sketch of the phase diagram is shown in FigEI With the 
help of Monte Carlo simulations we were able to bound 
the range of existence of the new phase. Unfortunately 
we have to admit that our results are not conclusive and, 
as discussed in the previous section, an alternative sce- 
nario is also possible. Nevertheless, we think that the 
existence of an intermediate phase is a very appealing 
possibility worth to being further investigated. 
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FIG. 17: A possible phase diagram of an array with Ts lattice. 
In the unfrustrated case (upper panel) we sketch the standard 
phase diagram which leads at T = to the SI transition. On 
the lower panel we present a possible scenario that emerges 
from our work. An new phase appears separating the normal 
from the superfluid phases. 



How is it possible to experimentally detect such a 
phase? In Josephson arrays, where one typically does 
transport measurement, the AB-insulator should be de- 
tected by looking at the temperature dependence of the 
linear resistance. On approaching the zero temperature 
limit, the resistance should grow as differently from 
the Mott insulating phase where it has an exponential ac- 
tivated behaviour. In optical lattices the different phases 
can be detected by looking at the different interference 
pattern (in the momentum density or in the fluctua- 
tionsS). A detailed analysis of the experimental probe 
will be performed in a subsequent publication. 

There are several issues that remain to be investigated. 
It would be important, for example, to see how the phase 
diagram of the frustrated system (and in particular the 
Aharonov-Bohm phase) is modified by a finite range of 
Uij and/or the presence of a finite chemical potential. 
An interesting possibility left untouched by this work is 
to study the fully frustrated array at uq = 1/2. In this 
case (for on-site interaction) the superfluid phase extends 



15 



down to vanishing small hopping. In this case a more 
extended AB insulating phase could be more favoured, 
and thus more clearly visible. 



Since the states \ip[Tk)) are eigenstates of Tit, the calcu- 
lation is reduced to the evaluation of the matrix elements 
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APPENDIX A: (2+l)D XY MAPPING 

We give here some of the technical details of the map- 
ping from the QPM to a (2 -f l)D - XY model. The 
latter one is particularly easy to be simulated numeri- 
cally: the state of the system and the effective action 
are both expressed in terms of phases on a ?>D lattice. 
Being n and canonically conjugated, it is possible to 
represent n as — and get the so-called quantum rotor 
Hamiltonian. For the sake of simplicity we consider a 
diagonal capacitance matrix. 

7i = Tit + T^u 

The partition function can be rewritten in a more con- 
venient way using the Trotter approximation: 

Z = Tr{(e-^(«'+«-))^^} 

= lim Tr{(e-'^^^^e-'^^"* -ho(AT2))^-}(A2) 

where Ht is imaginary time and At ~ [3 /L^ is the width 
of a time slice. The limit Ar must be taken to 

recover the underlying quantum problem. 

> 

Introducing complete sets of states \Lp{Tk)) with peri- 
odic boundary conditions on times (tq = = t^^ = 13) 
the trace can be written as 



Mrfe+i)|e'^-^|^(r,)). 



(A4) 



the matrix elements can be furtherly simplified going 
back to the charge representation (or angular momen- 
tum, since n is the generator of U{\) for the XY spin of 
a site): 



EH' 



'-[JJ? (,1- JU•^>i{Tk)-v^('rk+l)] ^ 



(A5) 



J" 



Using the Poisson summation formula, the sum over 
angular momentum configurations becomes a periodic se- 
quence of narrow gaussians around multiples of 27r 



n 



E 



27r 
C/At 



f [Vi(Tfc)-yi('rfe+l)-27rm]^ 



that is the Villain approximation to 



Tk 



n 



e V 



s['Pi{Tk)-Vi{Tk + l)] 



(A6) 



(A7) 



with dropped irrelevant prcfactors. 

What we get by means of this procedure is a map- 
ping of the QPM into a anisotropic classical (2 -I- 1)DXY 
model, with effective action S 



S 



[I - cos {(fi^k - 'Pj.k - 



{t.j).k 

i,{k, k') 

1 



tAr 



UAi 



(A8) 



(A9) 



where we used a symmetric notation for space and time 
lattice sites. Since critical properties are not expected to 
depend on the asymmetry of such model, and since for 
Ar — > we have Ksp — > 0, Kr — > oo with Ksp Kt = 
const., one can fix Ar = l/\/tU. It then follows that 
the coupling in the space and time directions are equal 
Ksp = Kt = K. The isotropic model, Ea. (|27|) is the one 
which is used in in the Monte Carlo simulations. 



fc=0 
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